Directional Ordering of Fluctuations in a Two-dimensional Compass Model 
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In the Mott insulating phase of the transition metal oxides, the effective orbital-orbital interaction 
is directional both in the orbital space and in the real space. We discuss a classical realization of 
directional coupling in two dimensions. Despite extensive degeneracy of the ground state, the 
model exhibits partial orbital ordering in the form of directional ordering of fluctuations at low 
temperatures stabilized by an entropy gap. Transition to the disordered phase is shown to be in the 
Ising universality class through exact mapping and multicanonical Monte Carlo simulations. 
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Recently, there has been growing interest in the effects 
of orbital degeneracy in the physics of transition metal 
oxide insulators 0, In these systems, the dominating 
energy scales for d-electrons on transition metal (TM) 
ions are the on-site Coulomb repulsion (which freezes 
out the charge degrees of freedom), the Hund's rule cou- 
pling, and the crystal field due to the surrounding oxygen 
ions. The latter two together determine the degenera- 
cies and degrees of freedom of spin and orbital on each 
transition metal ions. Spins and orbitals on neighbor- 
ing TM ions can then be coupled through the superex- 
change mechanism. In the case of orbitals, they are also 
coupled through the phonon-mediated cooperative Jahn- 
Teller mechanism 0. These couplings determine the low 
temperature properties of these systems. 

Because orbital coupling is intrinsically directionalQ, 
orbital ordering brings up some unusual questions. Es- 
pecially interesting is when the coupling along a given 
bond direction is Ising like, but with different Ising axes 
along different bond directions [4| . The Hamiltonian for 
orbitals is then given by 

H = -Jj2n ■ n^Tj .ny (1) 

W) 

where r is an isospin operator representing the orbital de- 
gree of freedom, and is a unit vector giving the Ising 
axis for the bond (ij) . For example, for e g orbitals in 
Perovskite structures, for the three different bond di- 
rections are coplanar and oriented relative to each other 
by 120°, giving rise to the so-called 120° Model0, H 0. 
This model is also applicable to t2 g orbitals on the three 
bonds of a honeycomb lattice as in planes of V2 03 0,0- 
On the other hand, for ti g orbitals on Perovskite struc- 
tures, the relevant model is the Compass Model [J UGH 
with Sjj — x, y, z for the three bond directions. The 
Compass Model may also feature as part of the spin- 
spin coupling of t2 g orbitals when spin-orbit interaction 
is taken into account II ll. A common feature of both the 



Compass Model and the 120° Model is the competition 
between bonds in different directions, with the result- 
ing frustration leading to macroscopic degeneracy of the 
classical ground state. 

In this Letter, we report analytical and numerical re- 
sults on a classical version of in 2D. The highly 
anisotropic coupling gives rise to interesting interplay be- 
tween one- and two-dimensional ordering, between con- 
tinuous and discrete spin physics, and between slow and 
fast modes. Our main result is that at low T > 0, there 
is no conventional ordering |l2j. but there is nevertheless 
LRO in the form of a directional ordering in fluctuations. 
This ordering corresponds to a partial breaking of the 4- 
fold symmetry and is stabilized by entropy. In this phase, 
the system exhibits spontaneous dimension reduction by 
essentially decoupling into one-dimensional (ID) chains 
running either horizontally or vertically along the bonds. 
Through exact mapping and extensive Monte Carlo sim- 
ulations, we show that this ordering transition belongs to 
the Ising universality class. We then discuss generaliza- 
tions of these results to the quantum case and in three 
dimensions, as well as their implications with respect to 
orbital ordering. 

Consider the classical Compass model on a square lat- 
tice of N = L x L sites, 

H = — J S y (SjxSj+x,X F SiySi + y^y), (2) 

i 

where Sj = (cosftj, sin^) represents either a real spin 
or an orbital isospin. On the square lattice, the sign of 
J can be gauged away, so we take J > 0. Along each 
row (column), we have a simple Ising model (IM) with 
quantization axis along x(y). Hamiltonian J2} has two 
types of discrete symmetry. I) There is a global four- 
fold rotation symmetry corresponding to simultaneously 
rotating the spins and lattice by multiples of 90°. II) 
In addition, it is also invariant under the ID spin flip 
transformation Si x — > —Si Xl Si y — > Si y for all i on any 
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one row and Si x — > S{ X , Si y — > — Si y for all i on any 
one column. Since symmetry I is two-dimensional (2D), 
we expect that it may be broken at finite T, while the 
fD nature of symmetry II should imply no symmetry 
breaking except possibly at T = 0. We will see indeed 
that this is the case, but the physics leading to it and 
their consequences are not trivially deduced from such 
symmetry considerations. 

The low temperature properties of J5J are further com- 
plicated by an additional 0(2) degeneracy of the ground 
state. Apart from a constant term, Eq. can be writ- 
ten as, 

H = — ^ [(cos0, - cos6> i+2 ) 2 + (sinfli - sm0 i+ y) 2 ] . 

i 

(3) 

Clearly 9i = is a ground state, as are the D = 2 x 2 L 
states obtained from it by the symmetry operations I 
and II. However, Eq. @ shows that the ground state 
energy is invariant under arbitrary global rotation of Qi = 
to 0i = 9. Unlike the isotropic XY model, where this 
invariance holds for each bond, here the energy loss from 
the horizontal bonds are compensated by energy gain in 
the vertical bonds. Thus, the ground state exhibits an 
0(2) degeneracy not related to the symmetries of H. Wc 
will see that this "accidental" degeneracy is lifted at finite 
temperatures by entropy due to slow mode fluctuations. 

Upon a redefinition of the spins through symmetry II, 
any of the ground states mentioned above can be re- 
cast as 8i = 9. To study slow mode physics, we start 
with the spin wave or harmonic approximation. Ex- 
panding (J3J to the second order in (pi = 6i — 9, we ob- 
tain the spin wave Hamiltonian in Fourier form, Hsw = 
JJ2k e k{9) |</5fc| 2 - For general 9, the spin wave spec- 
trum tk{9) = (1 — cos k x ) sin 2 9 + (1 — cosk y )cos 2 9 is 
an anisotropic 2D one, with zero modes at k x = k y = 0. 
However, for special ordering directions = 0, 7r/2, it, 
and 37r/2, the spectrum becomes ID like, independent of 
either k x (9 = 0, 7r) or k y (9 = tt/2, 3ir/2), and the den- 
sity of states of gapless excitations is ID rather than 2D. 
The high density of low lying states suggests an entropic 
mechanism to stabilize these four directions at T > 0. 

To put the discussion on a firmer footing, we employed 
a self-consistent harmonic approximation for the ordered 
phase at 9 = 0. Based on the Bogolyubov-Peierls theo- 
rem, F < Fq — (-Ho)o + (H)o, we compute the variational 
free energy using a trial Hamiltonian H Q = J^2 a k l^fel • 
Minimizing the free energy, we obtain 

a k = m + 7^,(1 - cosk x ) + j y (l - cosk y ) (4) 

where m is the self-consistent spin- wave gap, and j x and 
7j, are the self-consistent stiffness. 

At low temperatures, a ID spectrum with j x (T) = 0, 
7 y (T) = 1 - 0(T 2 / 3 ), and m(T) = ±T 2 / 3 + 0(T) is ob- 
tained. Within the self-consistent approach, anharmonic 
effects are incorporated into a shift of these parameters 



at finite T from their bare spin wave values. Most signifi- 
cantly we see that a gap m is generated, which suppresses 
the diverging ID fluctuations in the spin wave analysis, 
and stabilizes the ordering along one of the 4 special di- 
rections. 

To address whether there is ordering into one of the D 
degenerate ground states, we need to consider the effects 
of fast modes, or more precisely, abrupt spin flips. For 
this purpose, the continuous nature of the spins should 
not be crucial, so we discretize the Compass model into a 
"4-state Potts Compass model" with the same symmetry 
as (0), given by 

Hp = — J {Uju rtj-i-S t o- Oj CTj+x + n tT n i+y,T T i T i+y), (5) 

where on each site we have "occupation numbers" n ir7 = 
0, 1 and n iT = 1 — n ir7 . If n ia = 1, then there is an addi- 
tional internal degree of freedom a = ±1; and similarly 
for rii r and r. The correlation of these internal degrees of 
freedom with the occupation numbers together with the 
constraint in the latter couple these various variables. 

The partition function of the Potts compass mode 
takes the form Zp = Tr^^Tr/^ n \ exp(— f3Hp), where 
Tr' indicates that for a given configuration of {n icr } 
the trace over ct(t) should only be on those sites with 
rii a = 1(0). On the other hand, we note that if for ex- 
ample, riia — 0, then Hp is independent of <Ji, and trac- 
ing over Oi = ±1 simply gives a superfluous factor of 2. 
Thus, Tr' can be replaced by the unrestricted Tr to give 
Zp = 2- Ar Tr {n ^ } Tr {(Ti!ri} exp(-/3ff P ). The trace over a 
and r can now be easily done using transfer matrix since 
Hp consists of decoupled ID chains as far as a and r are 
concerned, resulting in Zp = Tr/ n4 a exp(— f3H e s), where 

H e ff = -Tln[cosh(/3J)] y^^(n i<7 n i+x ^ + n iT n i+ $ tT ) 

i 

-T ln[l + tanh L [fi J)] C a + ^ L> 7 ] , (6) 

a 7 

In the last two terms of H e s , C a — \\ i ni a ni +x ^, for all 
sites i in the row a, while D 1 — Y[j n jTnj+ ytT , for all sites 
j in the column 7. Under periodic boundary conditions, 
C a = 1 if rii a = 1 for all sites in row a and zero otherwise. 
Similarly, _D 7 = 1 if rii T = 1 for all sites in column 7 and 
zero otherwise. At T = 0, we recover the 2 x 2 L -fold 
degeneracy of the ground state due to symmetries I and 
II. 

At any T > 0, the last two terms in Eq. © are 
finite-sized terms that vanish in the thermodynamic limit 
L — > 00. Ignoring them, we may rewrite H e g in terms of 
n ia ,r = 3(1 T Hi) as 

H cS = -2NJ - jy^(^/i t+ £ + Hi[i t+y ). (7) 

i 

The 4-state Potts Compass model is thus mapped exactly 
into the 2D Ising Model (2DIM). The coupling constants 
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of the two models are related by J = Tln[cosh(J/T)]/4. 
From the 2DIM exact f c = 2J/ln(l + \/2), we conclude 
that the Potts Compass model has long-ranged order 
(LRO) for all T < T c = 0.4084J. What is the nature 
of this LRO? First, note that because the trace over a 
and t are for decoupled chains, (cr,) and (n) = for all 
T > 0. Instead, the 2DIM uniform ordering of (/Xj) corre- 
sponds to (n a ) — (n T ) ^ 0. In other words, the ordering is 
not of the spins spontaneously pointing along one of the 4 
possible states, but of them having stronger fluctuations 
in 2 of the 4 states, henceforth called directional order- 
ing of fluctuations. In this phase, the Z4 symmetry of 
the Compass model is only partially broken into Z2 x ^2- 
While at T = 0, the ground state has macroscopic degen- 
eracy, the free energy has only two degenerate minima at 
T = 0+ 

Based on symmetry considerations and on our entropy 
stabilization arguments earlier, we expect the above con- 
clusions to hold also for the continuous Compass model 
except for the value of T c . To confirm this and to rule 
out a preemptive first order transition, we perform Monte 
Carlo simulations. However, such simulations are compli- 
cated by the strong size dependence that originates from 
the finite-sized terms in Eq. ijfjj under periodic bound- 
ary condition when the ID correlation length £id exceeds 
the linear system size L at sufficiently low temperature. 
This boundary contribution, however, can be eliminated 
if, for each row and column, we pick one bond and allow 
the sign of the coupling constant on that bond to be an 
additional degree of freedom. The data presented below 
are obtained under this "annealed" boundary condition. 
In view of the large number of degenerate low energy 
states arising from Symmetry II, the simulation is per- 
formed using the multicanonical Monte Carlo sampling 
scheme |l5|. Details will be published elsewhere. 

Figure 1(a) shows the average of the directional 
bond order parameter q = N -1 5^(5?, ~ &iy) = 
N- 1 J2t cos 29 t against T for L = 8, 12, 16, 24, 32 and 48. 
Existence of an ordered state at low temperatures is ev- 
ident from the data. To locate the transition point, we 
computed the Binder cumulant B = 1 — (q 4 ) / (3(q 2 ) 2 ) for 
various system sizes, as shown in Fig. 1(b). From the 
crossing of the curves, we estimate T c = 0. 147 ± 0.001 for 
the continuous Compass model. The value of B at the 
crossing point for large system sizes also agrees reason- 
ably well with the 2DIM result B c = 0.61069 . . ..[iH 

Figure 2 shows the specific heat data for the six differ- 
ent sizes up to L — 48. A weak divergence near the T c 
value determined above is clearly seen. The inset shows 
the same data in the critical region after subtraction of a 
background linear function, plotted using the scaled vari- 
ables according to the finite-size scaling form derived in 
Ref.[I3. Apart fr om the smallest size at L — 8, the data 
collapse is quite satisfactory. 

The simulation data shows quite convincingly that di- 
rectional bond ordering exists at low temperatures, and 
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FIG. 1: (a) Directional bond ordering parameter against tem- 
perature T for the 2D Compass model at various system sizes, 
(b) Binder cumulant against temperature. The critical tem- 
perature is estimated at T c = 0.147 ± 0.001. 




FIG. 2: Specific heat against temperature for the 2D Compass 
model under the annealed boundary condition. Inset: same 
data plotted using scaled variables. 



the transition to the disordered phase is in the Ising uni- 
versality class. The transition temperature of the contin- 
uous Compass model, on the other hand, is considerably 
lower than that of the Potts Compass model. We at- 
tribute this to the softening of domain wall energy in the 
continuous model. Indeed, the free energy gap between 
the favored orientation and other states does not remain 
constant, but in fact vanishing as T 2 / 3 at low T. 

We next discuss the implications of our results to or- 
bital ordering, taking the spin in the Hamiltonian to rep- 
resent the orbital isospin. For that we should take into 
account these spins will then have S < 00, i.e. quan- 
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turn mechanical rather than classical, and so there will 
be fluctuations about the classical ground state even at 
T = 0. For TMO, the relevant S are 1/2 and 1 for 
double and triple orbital degeneracy respectively. For 
slow modes, the physics is somewhat similar to that of 
the Heisenberg antiferromagnet on the Kagome lattice, 
and also similar to straightforward extension of the finite 
temperature classical analysis above to 2 + 1 dimensions, 
with renormalized spin wave fluctuations lifting the 0(2) 
degeneracy while also generating a gap that stabilizes 
the 9 — 0, 7r/2, 7r, and 3n/2 ordered states. The same 
2 + 1 dimension argument would imply that the fast mode 
physics would not destroy the conventional ordering at 
T = 0. However, the true description of the fast mode 
physics would involve understanding the physics of in- 
stanton tunneling of the quantum Compass Model, and 
is beyond the scope of this paper. At any rate, at low 
but finite T, we expect both the classical and quantum 
mechanical models to exhibit directional but not conven- 
tional ordering, thereby leaving the orbital degeneracy 
completely (S = 1/2) or partially (S = 1) unbroken. 
In the case of S — 1, the directional ordering reduces 
the triple orbital degeneracy to double degeneracy. In 
the case of S — 1/2, Sf x = Sf = 1/4, and the direc- 
tional order parameter q defined above = 0. Instead, we 
define an alternative directional order parameter appli- 
cable for all S, namely, r =< S ix S i+XiX - S iy Si + y^ >, 
which like q shows LRO in the direction of fluctuations in 
isospin space. An advantage of r over q is that it explic- 
itly displays the energy difference between horizontal and 
vertical bonds and hence the broken lattice rotation sym- 
metry. A consequence of this is that, when the couplings 
of the orbital isospin to lattice modes are included, the 
directional ordering will be necessarily accompanied by 
a lattice distortion so that the bond length in horizontal 
and vertical directions become unequal. 

Our results can be generalized to the Compass model 
on the 3D cubic lattice, with now a 3-component spin 
S = (S x , S y , S z ). The slow mode physics is basically 
identical to the 2D case, so that orderings along ±x, 
±y, ±z have ID spin wave spectrum within harmonic 
approximation. Anharmonic terms then generate an en- 
tropy gap as in the 2D case. Restricting to the 6 spe- 
cial directions, we can then consider the corresponding 
6-state Potts compass model 18] . The mapping we used 
to map the 2D 4-state Potts Compass model into the 
2D Ising model can be applied here too. However, in 
this case, the mapping does not result in an ordinary 
3-state Potts model, but in a 3-state Potts Compass 
model H = niani +x ^+ni T ni + y :T +n itl n i+x ^ 7 where 
f^ia,r,n — 0, 1 and rii a + 7ii T + rii^ = 1. There is no known 
solution to this model. Heuristic domain wall analysis 
suggests that directional ordering is stable at low T, but 
more rigorous calculation is necessary for this to be con- 
clusive. Unlike the 2D case, the degeneracy of the di- 
rectional ordered state, while less than the ground state 



degeneracy of D = 3 x (2 x 2 L ) L , is also macroscopically 
degenerate. For example, beginning with the directional 
ordered state rii a = 1, one can take all the spins on an ar- 
bitrary number of planes perpendicular to the z-axis and 
convert them all to 7ij T = 1 with no cost in free energy. 
The degeneracy is then found to be 3 x 2 x 2 L . Notice 
that the fluctuation effects are stronger in 3D than 2D, 
and in fact, in general, the higher the dimension, the 
stronger the fluctuations, a point previously pointed out 
by Khomskii and Mostovoy in Ref . 0] . 

In conclusion, we have established that the Compass 
model has a low temperature phase characterized by 
(Si) = 0. but with long-ranged correlations in the di- 
rection of fluctuations in both isospin and lattice space. 
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